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I. INTRODUCTION 



Various methods have been proposed to determine the 
Unitarity Triangle angle 7 [U, 0, Q of the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix Q us- 
ing i)i*)0]^{*)- decays, where the symbol D^*^^ in- 
dicates either a 1?*^*^° or a meson. A B can decay 
into a D^*^'^ K^*^~ final state via a6— >cora6— s-w medi- 
ated process and CP violation can be detected when the 
D^*)" and the _D(*)o decay to the same final state. These 
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processes are thus sensitive to 7 = a.i:g{—V*i^Vud/V*^Vcd}- 
The present determination of 7 comes from the combina- 
tion of several results obtained with the different meth- 
ods. In particular, the Dalitz technique when used 
to analyze B' I)(*)Oif(*)- decays, is very powerful, 
resulting in an error on 7 of about 24° and 13° for the 
BaBAR and Belle analyses respectively ( P, 'g'] ) . These re- 
sults are obtained from the simultaneous exploitation of 
the three decays of the charged B mesons {B^ D'^K^, 
D*°K- and D°K*-) and, in the case of BaBAR, from 
the study of two final states for the neutral D mesons 
{Ksn+T:~ and KsK+R-). 

In this paper we present the first measurement of the 
angle 7 using neutral B meson decays. We reconstruct 
S° b°K*°, with K*° K+n- (charge conjugate 
processes are assumed throughout the paper and K*^ 
refers to if* (892)°), where the flavor of the B meson is 
identified by the kaon electric charge. Neutral D mesons 
are reconstructed in the KsTT'^Tr~ decay mode and are 
analyzed with the Dalitz technique The final states 
we reconstruct can be reached through 5 — > c and 6 — > u 
processes with the diagrams shown in Fig. [TJ The corre- 
lation within the flavor of the neutral D meson and the 
charge of the kaon in the final state allows for discrim- 
inating between events arising from 6 — > c and 6 — s- u 
transitions. In particular it is useful for the following 
discussion to stress that 6 — > u (5°) transitions lead to 



D^K'^TT final states and h 
D°K-TT+ final states. 



/-ftO 



u {B ) transitions lead to 



corresponding to the if *° . For this purpose we introduce 
the quantities rs, k, and Ss defined as 



_ riB^^D°K+n-) _ JdpAl{p) 



r2 = 

I Q — 



ke 



_ J dp AMMp)e''^''^ ^ 
^JdpA^M JdpAlip)' 



(2) 



where < fc < 1 and 63 € [0, 27r]. The amplitudes for 
the b c and b u transitions, Adp) and A^ip), are 
real and positive and d{p) is the relative strong phase. 
The variable p indicates the position in the D^K'^tt~ 
Dalitz plot. In case of a two-body B decay, rs and 5s 
become rs — \Au\/\Ac\ and 5b (the strong phase differ- 
ence between A^ and A^) and k = 1. Because of CKM 
factors and the fact that both diagrams, for the neutral 
B decays we consider, are color-suppressed, the average 
amplitude ratio rs in B^ D^K*^ is expected to be in 
the range [0.3, 0.5], larger than the analogous ratio for 
charged B^ D^K^ decays (which is of the order of 
10% [1, Q). An earlier measurement sets an upper hmit 
rs < 0.4 at 90% probability [l^]- A phenomenological ap- 
proach [ll| proposed to evaluate rs in the B^ D^K^ 
system gives =0.27±0.18. 



II. EVENT RECONSTRUCTION AND 
SELECTION 




FIG. 1: Feynman diagrams for the decays B'^ D'^K*^ (up 



left, b ^ c transition), — > D" K*" (up right, b ~* u tran- 
jjOj^tO (•jJq^jj igf^^ ^ _^ f. transition), and 

u transition). A K*" is a 
decay product of a while a K*" results from a B" decay. 



sition), b'^ 

3° D°K*° (down right, b 



When analyzing — > D^K*^ decays, the natural 
width of the K*^ (50 MeV/c^) has to be considered. In 
the K*^ mass region, amplitudes for decays to higher- 
mass Kt: resonances interfere with the signal decay am- 
plitude and with each other. For this analysis we use 
effective variables, introduced in Ref. 0], obtained by in- 
tegrating over a region of the B^ — > D^K^'k~ Dalitz plot 



The analysis presented in this paper uses a data sample 
of 371 X 10*^ BB pairs collected with the BaBAR detec- 
tor at the PEP-II storage ring. Approximately 10% of 
the collected data (35 fb^-'^) have a center-of-mass (CM) 
energy 40 MeV below the T{AS) resonance. These "off- 
resonance" data are used to study backgrounds from con- 
tinuum events, e+e^ ~^ QQ il ~ d, s, or c). 

The BaBAR detector is described elsewhere [l^ . 
Charged-particle tracking is provided by a five-layer sil- 
icon vertex tracker (SVT) and a 40-layer drift cham- 
ber (DCH). In addition to providing precise position 
information for tracking, the SVT and DCH also mea- 
sure the specific ionization (dE/dx), which is used for 
particle identification of low-momentum charged parti- 
cles. At higher momenta {p > 0.7 GeV/c) pions and 
kaons are identified by Cherenkov radiation detected in 
a ring-imaging device (DIRC) . The position and energy of 
photons are measured with an electromagnetic calorime- 
ter (EMC) consisting of 6580 thallium-doped Csl crys- 
tals. These systems are mounted inside a 1.5T solenoidal 
super-conducting magnet. 

We reconstruct B° D°K*^ events with K*° ^ 
K^TT^ and i^sTr+vr^. The event selection, de- 

scribed below, is developed from studies of off-resonance 
data and events simulated with Monte Carlo techniques 
(MC). The Ks is reconstructed from pairs of oppositely- 
charged pions with invariant mass within 7 MeV/c^ of 
the nominal Ks mass 13|, corresponding to 2.8 standard 
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deviations of the mass distribution for signal events. We 
also require that cosaKsiD^) > 0.997, where a^f 
is the angle between the Kg line of flight (line between 
the and the Kg decay points) and the Ks momen- 
tum (measured from the two pion momenta) . Neutral D 
candidates are selected by combining Kg candidates with 
two oppositely-charged pion candidates and requiring the 
D'^ invariant mass to be within 11 MeV/c^ of its nomi- 
nal mass [l^, corresponding to 1.8 standard deviations 
of the mass distribution for signal events. The Kg and 
the two pions used to reconstruct the D*^ are constrained 
to originate from a common vertex. The charged kaon is 
required to satisfy kaon identification criteria, which are 
based on Cherenkov angle and dE/dx measurements and 
are typically 85% efficient, depending on momentum and 
polar angle. Misidentification rates are at the 2% level. 
The tracks used to reconstruct the K*'^ are constrained 
to originate from a common vertex and their invariant 
mass is required to lie within 48 MeV/c^ of the nominal 
K*" mass [l3|. We define OhcI as the angle between the 
direction of fiight of the charged K in the K*^ rest frame 
with respect to the direction of fiight of the K*^ in the 
B rest frame. The distribution of cos 9 Hei is expected to 
be proportional to cos^ Ouei for signal events, due to an- 
gular momentum conservation, and flat for background 
events. We require \ cos9Hei\ > 0.3. The cuts on the 
if*" mass and on | cos Onei \ have been optimized maxi- 
mizing the function S'/V S + B, where S and B are the 
expected numbers of signal and background events re- 
spectively, based on MC studies. The 5° candidates are 
reconstructed by combining one 13° and one if*" can- 
didate, constraining them to originate from a common 
vertex with a probability greater than 0.001. The dis- 
tribution of the cosine of the B polar angle with respect 
to the beam axis in the e'^e" CM frame, cosOb, is ex- 
pected to be proportional to 1 — cos^ 9b- We require 
I cos 6*5 1 < 0.9. 

We measure two almost independent kinematic vari- 
a bles: the beam-energy sub stituted mass toes = 
-phY/El -pB^, and the energy differ- 
ence AE = Eg — Eq/2, where E and p are energy and 
momentum, the subscripts B and refer to the candidate 
B and to the e+e~ system respectively and the asterisk 
denotes the e+e~ CM frame. For signal events, TOes is 
centered around the B mass with a resolution of about 2.5 
MeV/c^ , and AE is centered at zero with a resolution of 
12.5 MeV. The B candidates are required to have AE in 
the range [—0.025, 0.025] GeV. As it will be explained in 
Sec. IIV| the variable toes is used in the fit procedure for 
the signal extraction. For this reason, the requirement 
on it are quite loose: toes € [5.20,5.29] GeV/c^. The 
region 5.20GeV/c^ < toes < 5.27 GeV/ c^, free from any 
signal contribution, is exploited in the fit to characterize 
the background directly on data. The proper time inter- 
val At between the two B decays is calculated from the 
measured separation, Az, between the decay points of 
the reconstructed B {B^^c) and the other B {Both) along 
the beam direction. We accept events with calculated At 



uncertainty less than 2.5 ps and |At| <20 ps. In less than 
1% of the cases, multiple candidates are present in the 
same event and we choose the one with reconstructed D° 
mass closest to the nominal mass [l^ . In the case of two 
B candidates reconstructed from the same D^, we choose 
the candidate with the largest value of | cos 9Hei \ ■ The 
overall reconstruction and selection efficiency for signal, 
evaluated on MC, is (10.8 ± 0.5)%. 



III. BACKGROUND CHARACTERIZATION 

After applying the selection criteria described above, 
the background is composed of continuum events 
(e+e^ qq, q — u, d, s, c) and T(4S) — > BB events 
( "BB" , in the following) . To discriminate against the 
continuum background events (the dominant background 
component), which, in contrast to BB events, have a jet- 



like shape, we use a Fisher discriminant T [ij]. The 
discriminant J-' is a linear combination of three vari- 
ables: cos 9 thrust , the cosine of the angle between the 
B thrust axis and the thrust axis of the rest of the event, 
= J2iPi' and L2 = X^iPil cos^il^. Here, pi is the 
momentum and 9i is the angle with respect to the thrust 
axis of the B candidate. The index i runs over all the re- 
constructed tracks and energy deposits in the calorimeter 
not associated with a track. The tracks and energy de- 
posits used to reconstruct the B are excluded from these 
sums. All these variables are calculated in the e+e" CM 
frame. The coefficients of the Fisher discriminant, chosen 
to maximize the separation between signal and contin- 
uum background, are determined using signal MC events 
and off-resonance data. A cut on this variable with 85% 
efficiency on simulated signal events would reject about 
80% of continuum background events, as estimated on 
off-resonance data. We choose not to cut on the Fisher 
discriminant, as we will use this variable in the fit pro- 
cedure to extract the signal. The variable At gives fur- 
ther discrimination between signal and continuum events. 
For events in which the B meson has been correctly re- 
constructed, the At distribution is the convolution of a 
decreasing exponential function e^*/"^^ (with tb equal to 
the B lifetime) with the resolution on Az from the detec- 
tor reconstruction. The distribution is then wider than in 
the case of continuum events, in which just the resolution 
effect is observed. 

The i^ioc decay point is the common vertex of the B 
decay products. The i^oth decay point is obtained using 
tracks which do not belong to B^^c and imposing con- 
straints from the i^icc momentum and the beam-spot lo- 
cation. 

Background events for which the reconstructed Ks, 
TT"*", and TT^ come from a real Z)" ("true i)"", in the 
following) are treated separately because of their distri- 
bution over the D° Dalitz plane. A fit to the Kstt~^7t~ 
invariant mass distribution for events in the toes side- 
band (toes < 5.27 GeV/ c^) has been performed on data 
to obtain the fraction of true D" equal to 0.289 ± 0.028. 
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This value is in agreement with that determined from 
simulated BB and continuum background samples. 

Background events with final states containing 
D^h^n~ or D^h'TT^, where is a candidate 
and Kstt^tt^, can mimic b ^ u mediated sig- 

nal events (see Fig. [Ij. The fraction of these events 
(relative to the number of true events), defined as 

nb^u — j^(^o0h + TT-)+N{D0h-7T+)+N{D'>h+Tr-} + N(D0h-TT+y 

has been found to be 0.88 ± 0.02 and 0.45 ± 0.12 in BB 
and continuum MC events respectively. 

Studies have been performed on B decays, which have 
the same final state reconstructed particles as the signal 
decay (so called peaking background). From MC stud- 
ies, we identify three possible background sources of this 
kind: B° D°K*° {K*° ^ K+tt-, D° 7r+7r-7r+7r-), 
B° D°p° (pO _> TT+TT-, D° KsTT+TT-, where p° is 
reconstructed as a K*'^ with a misidentified pion) and 
charmless events of the kind 5° K*°KsKs- To 
precisely evaluate the selection efficiency for and 
b°K*° with b° 7r+7r~7r+7r-, dedicated MC sam- 
ples have been generated, resulting in (0.04 ± 0.02)% or 
(0.18 ± 0.04)%, respectively. With these efhciencies, we 
expect to select about 0.9 p^ events and 0.1 Att 
events in 371 x 10^ BB pairs. In the latter case the re- 
quirement on aKs rejects most of the background, while 
for p^ the cuts on A_E and the particle identifica- 
tion of the are the most effective. The number of 
charmless background events has been evaluated on data 
from the D'^ mass sidebands, namely Mj^g^+^- in the 
range [1.810, 1.839] or [1.889, 1.920] GeV/c^; we obtain 
Npeak = —5 ± 7 events, consistent with 0. Hence we as- 
sume these background sources can be neglected in our 
signal extraction procedure; the effects of this assump- 
tion are taken into account in the evaluation of the sys- 
tematic uncertainties. The remaining BB background is 
combinatorial. 



IV. LIKELIHOOD FIT AND MEASURED YIELD 

We perform an unbinned extended maximum likeli- 
hood fit to the variables tties, ^ and At, in order to 
extract the signal, continuum and BB background yields, 
probability density function (PDF) shape parameters, 
and CP parameters. We write the likelihood as 

a i—1 

where P"(i) and Na are the PDF for event i and the 
total number of events for component a (signal, BB 
background, continuum background). Here N is the 
total number of selected events and 77 is the expected 
value for the total number of events, according to Pois- 
son statistics. The PDF is the product of a "yield" 
PDF V°'{mE.s)V°'{T)V"'{M) (written as a product of ID 
PDFs since mES, ^ and At are not correlated) and of 



the Dalitz plot dependent part: V"{rn?^,m'^_) (where 
to5, = TO?, 4. and m? — to?, _ ) . 

The Toes distribution is parametrized by a Gaussian 
function for the signal and by an Argus function Ts'] that 
is different for continuum and BB backgrounds. The T 
distribution is parametrized using a asymmetric Gaus- 
sian distribution for the signal and BB background and 
the sum of two Gaussian distributions for the continuum 
background. For the signal, |At| is parametrized with an 
exponential decay PDF e~*/'^ in which r = t^o [l^, con- 
volved with a resolution function that is a sum of three 
Gaussians |T6j . A similar parametrization is used for the 
backgrounds using exponential distributions with effec- 
tive lifetimes. 

The continuum background parameters are obtained 
from off-resonance data, while the BB parameters are 
taken from MC. The fractions of true I?" and the ratios 
Rb^u in the backgrounds are fixed in the fit to the values 
obtained on data and MC respectively. 

Using the effective parameters defined in Eq. [1] and 
m the partial decay rate for events with true 13° can be 
written as follows: 

r(S° ^ D[Ksn-TT+]K+TT-) cx \V^J^\^ + 
rl\r^;'\^ + 2krs\V^_:^\\Vf'\ C08{Ss + 6d- - 7) , (4) 
T{W ^ D[Ksn-7r+]K-T:+) cx \r+^\^ + 
rl\V^_:^\^ + 2fcrs|7'f IIT'^'^I cos(<55 + Sd+ + 7) , (5) 

where V^^ = V^'^ (ml , mi) , V^^ = V^'^imi,ml) and 
where Sd+ = '^d(™+, to^) is the strong phase difference 
between 7^^'^ and and Sd- = <5d(to^,to^) is the 

strong phase difference between and V^^. 

For the resonance structure of the ~> ifgTr+Tr^ 
decay amplitude, V^^, we use the same model as docu- 
mented in |5|] . This is determined on a large data sample 
(about 487000 events, with 97.7% purity) from a Dahtz 
plot analysis of mesons from D*'^ D'^tt~^ decays 
produced in e+e" cc events. The decay amplitude is 
parametrized, using an isobar model, with the sum of the 
contributions of ten two-body decay modes with interme- 
diate resonances. In addition, the K-matrix approach [iTj 
is used to describe the S-wave component of the tt+tt" 
system, which is characterized by the overlap of broad 
resonances. The systematic effects of the assumptions 
made on the model used to describe the decay ampli- 
tude of neutral D mesons into Kstt^tt~ final states are 
evaluated, as it will be described. To account for possi- 
ble selection efficiency variations across the Dalitz plane, 
the efficiency is parametrized with a polynomial function 
whose parameters are evaluated on MC. This function is 
convoluted with the Dalitz distribution V^^. The distri- 
bution over the Dalitz plot for events with no true D'^ is 
parametrized with a polynomial function whose parame- 
ters are evaluated on MC. 

Following Ref. [11] , we have performed a study to eval- 
uate the possible variations of rg and k over the 
D^K^'K~ Dalitz plot. For this purpose we have built 
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a -B Dalitz model suggested by recent measurements 
[H m, including K*{892f, K*{U30)°, K*{U30)°, 
iC*(1680)0, L>^2(2573)±, £'^(2460)*, and d{{230S)^ 
contributions. We have considered the region within 48 
MeV/c^ of the nominal mass of the K*{892)'^ resonance 
and obtained the distribution of rs and k by randomly 
varying all the strong phases ([0,27r]) and the amplitudes 
(within [0.7,1.3] of their nominal value). The ratio be- 
tween 6 — > M and 6 — s- c amplitudes for each resonance 
has been fixed to 0.4. In the K*{892)^ mass region, we 
find that rg varies between 0.30 and 0.45 depending upon 
the values of the contributing phases and of the ampli- 
tudes. The distribution of k is quite narrow, centered at 
0.95 with an r.m.s. of 0.02. The study has been repeated 
varying the ratio between b u and b c amplitudes 
between 0.2 and 0.6, leading to very similar results. For 
these reasons the value of k has been fixed to 0.95 and 
a variation of 0.03 has been considered for the system- 
atic uncertainties evaluation. On the contrary, rg will be 
extracted from data. 

We perform the fit for the yields on data extracting 
the number of events for signal, continuum and BB, as 
well as the slope of the Argus function for the BB back- 
ground. The fitting procedure has been validated using 
simulated events. We find no bias on the number of fitted 
events for any of the components. The fit projection for 
WES is shown in Fig. H We find 39±9 signal, 231±28 BB 
and 1772 ± 48 continuum events. In Fig. [5] we also show, 
for illustration purposes, the fit projection for tties, after 
a cut on !F > 0.4 is applied, to visually enhance the sig- 
nal. Such a cut has an approximate efficiency of 75% on 
signal, while it rejects 90% of the continuum background. 



V. DETERMINATION OF 7 



used to describe true D" Kstt~^tt~ decays is evalu- 
ated on data by repeating the fit with models alternative 
to the nominal one, as described in detail in The 
D'^ KsTT^TT~ Dalitz model is known to be the source 
of the largest systematic contribution in this kind of mea- 
surements [q . All the other contributions have been 
evaluated on a high statistics simulated sample in order 
not to include statistical effects. To evaluate the con- 
tribution related to mESi ^ and At PDFs, we repeat 
the fit by varying the PDF parameters obtained from 
MC within their statistical errors. To evaluate the un- 
certainty arising from the assumption of negligible peak- 
ing background contributions, the true 13° fraction and 
Rh~tu in the background, we repeat the fit by varying the 
number of these events and fractions within their statis- 
tical errors. The uncertainty from the assumptions on 
the factor k is also evaluated. The reconstruction effi- 
ciency across the Dalitz plane for true events and the 
Dalitz plot distributions for background with no true 
have been parametrized on MC using polynomial func- 
tions. Systematics uncertainties have been evaluated by 
repeating the fit assuming the efficiency and the distri- 
bution for these backgrounds to be flat across the Dalitz 
plane. 



Systematics source 


A7[°] 




Ars(10-^) 


Dalitz model for signal 


6.50 


15.80 


6.00 


PDF shapes 


1.50 


2.50 


5.20 


Peaking background 


0.14 


0.12 


0.04 


k parameter 


0.07 


1.20 


7.10 


True D" in the background 


0.05 


0.03 


1.00 


Rb-tu 


0.01 


1.10 


1.90 


Efficiency variation 


0.31 


0.62 


0.61 


Dalitz background param. 


0.03 


0.27 


0.20 


Total 


6.70 


16.10 
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From the fit to the data we obtain a three-dimensional 
likelihood C for 7, ^5 and rs which includes only sta- 
tistical uncertainties. We convolve this likelihood with a 
three-dimensional Gaussian that takes into account the 
systematic effects, described later, in order to obtain the 
experimental three-dimensional likelihood for 7, 5s and 
rs- From simulation studies we observe that, due to the 
small signal statistics and the high background level, rs is 
overestimated and the error on 7 is underestimated, when 
we project the experimental three-dimensional likelihood 
on either rs or 7, after integrating over the other two 
variables. This problem disappears if either rs is fixed in 
the fit or if we combine the three-dimensional likelihood 
function (7,(5s,''s) obtained from this data sample with 
external information on rg. In the following we will show 
the results of both these approaches. 

The systematic uncertainties, summarized in Table HI 
are evaluated separately on 7, (5s and rs and consid- 
ered uncorrelated and Gaussian. It can be noted that 
the systematic error is much smaller than the statistical 
one. The systematic uncertainty from the Dalitz model 



TABLE I: Systematics uncertainties on 7 , (5s , and rs- 

In Fig. [31 we show the 68% probability region obtained 
for 7 assuming different fixed values of rs and integrat- 
ing over ^5. For values of rg < 0.2 we do not have a 
significant measurement of 7. The value of (the fixed) rg 
does not affect the central value of 7, but its error. For 
example, for rs fixed to 0.3, we obtain 7 = (162 ± 51)°. 
On MC, for the same fit configuration, the average er- 
ror is 45° with a r.m.s. of 14°. The BABMi analysis for 
charged B decays Hj, using the same Dalitz technique 
for D'^ — !■ A'sTT+TT^, gives, for a similar luminosity, an er- 
ror on 7 of 29°, from the combination of B^ D^K^, 
B± D*°K^ and B± ^ D^K*^. The use of neutral B 
decays can hence give a contribution to the improvement 
of the precision on 7 determination comparable with that 
of a single charged B channel. 

Combining the final three-dimensional PDF with the 
PDF for rs measured with an ADS method 0], recon- 
structing the neutral D mesons into flavor modes [loj . 
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FIG. 2: rriES projection from the fit (a). The data are indicated with dots and error bars and the different fit components are 
shown: signal (dashed), BB (dotted) and continuum (dot-dashed). With a different binning (b), mss projection after a cut on 
JT > 0.4 is applied, to visually enhance the signal. and At projection from the fit (c, d). 
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while, at 95% probability: 
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FIG. 3: 68% probability regions obtained for 7, for different 
values of rs- For values of rs lower than 0.2, the distribution 
obtained for 7 is almost flat and hence does not allow to 
determine significative 68% probability regions. The solution 
corresponding to a 180° ambiguity is not shown. 



7 e [77,247]° or [257,426]°; 
6s € [-23,147]° or [157,327]°; 
rs < 0.55. 



(9) 
(10) 
(11) 



The preferred value for 7 is somewhat far from the 
value obtained using charged B decays, which is around 
75° for both BaBAR and Belle Dalitz analyses, but is 
compatible with both the results within about 1.5f7. In 
Fig. m we show the distributions we obtain for 7, rg 
and 7 vs. rs (the 68% and 95% probability regions are 
shown in dark and light shading respectively). The one- 
dimensional distribution for a single variable is obtained 
from the three-dimensional PDF by projecting out the 
variable and integrating over the others. 



we obtain, at 68% probability. 



VI. CONCLUSIONS 



In summary, we have presented a novel technique for 

/^ CO _L c;fi^o foAoj- Ka\° /a\ extracting the angle 7 of the Unitarity Triangle in B° 

7 = (162 ± 56) or (342 ± 56) ; (6) ^0^*0 (^go ^ f)'>K*°)with the ^ K+tt- (/?*« ^ 

Ss = (62 ±57)° or (242 ±57)°; (7) K-Tr+), using a Dalitz analysis of i)o ^ KsTt+tt-. With 

rs < 0.30; (8) the present data sample, interesting results on 7 (Eqs. 
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FIG. 4: 68 % (dark shaded zone) and 95 % (light shaded zone) probabihty regions for the combined PDF projection on the 7 
vs rs plane (a), 7 (b) and rs (c). 



[6l [9]) and rs (Eqs. [H [TT|l are obtained when combined 
with the determination of rs from the study of decays 
into flavor modes. The result for 7 is consistent, within 
1.5(7, with the determination obtained using charged B 
mesons. If the ratio rs is found to be of the order of 0.3, 
the use of neutral B mesons, proposed here, could give 
results on 7 as precise as those obtained using similar 
techniques and charged B mesons 
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